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Abstract 

We argue that there is a new hquid phase in the two-dimensional electron 
system in Si MOSFETs at low enough electron densities. The recently ob- 
served metal-insulator transition results as a crossover from the percolation 
transition of the liquid phase through the disorder landscape in the system 
below the liquid-gas critical temperature. The consequences of our theory are 

discussed for variety of physical properties relevant to the recent experiments. 
PACS numbers: 71.30.-hh,73.40.Qv,73.40.Hm 
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The scaling theory of locahzation [|l| of non-interacting electrons tells us that disorder 
is a relevant perturbation in two-dimension: The system is an insulator at low enough tem- 
peratures. The recent discovery of a possible metal-insulator transition in Si MOSFETs by 
Kravchenko et al. apparently defies this long-held belief. This should not be surprising be- 
cause the dominant Coulomb interaction in Si MOSFETs may invalidate the non-interacting 
scaling theory. In a typical Si MOSFET sample where the transition is observed, the ratio 
between the average Coulomb interaction and the Fermi energy is about 20. These intrigu- 
ing experiments generate renewed interests in the effects of interaction in disordered 
systems One of the striking features of the system is the I-V nonlinearity at average 
electric fields much weaker than the expected value determined by the effective temperature 
of the electrons . This raises questions of whether the system can be described by a theory 
for homogeneous systems. (See the discussion on electric field dependence in the text for 
more detail.) As we argue in this letter, the mobile positive background in Si MOSFETs 
allows macroscopic inhomogeneity to occur more easily at low enough average electron den- 
sities. The metal- insulator transition, the small field nonlinear I-V, as well as a host of other 
phenomena result from the combined effect of the inhomogeneity and the disorder in the 
system. 

In a remote doped or modulation doped GaAs/AlGaAs sample, the electron system is 
well described by the classic jellium model of electron gas. The roughly uniform positive 
background is fixed by the bulk crystal and it does not participate in the dynamics of the 
electron system. On the other hand, in an inverted Si MOSFET, the positive background 
charges at the metal-oxide interface are mobile. They couple to the dynamics of the electron 
system, as exhibited by the recent effective mass measurements in such systems [0. This 
allows for more possibilities for the physical properties of the electron system in Si MOS- 
FETs. However, we want to emphasize that the mobile positive background is not absolutely 
essential for the macroscopic inhomogeneity to occur. A slowly-varying disorder potential 
combined with electron-electron correlation, for example, may also favor macroscopic inho- 
mogeneity in the system. 
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We now consider what happens at low enough electron densities in a Si MOSFET. For 
the convenience of our argument, we treat the mobile positive background charges as holes. 
We consider the adiabatic evolution of the system as the oxide thickness increases. At zero 
thickness when the electrons and holes are at the same plane, the phase diagram of such 
an electron- hole system, shown in Fig. is well established 1^,0]. The high density phase 
is a metallic electron-hole liquid. The low density phase is an exciton gas. There is a two- 
phase coexistence region below the critical temperature Tc. As the oxide thickness increases 
beyond the Bohr radius as, the excitons in the gas phase dissolve into separated electrons 
and holes. In a disordered sample, these electrons occupy the lowest lying localized states in 
the system. Because these states are strongly localized, the internal degrees of freedom, such 
as spin and valley indices, do not play important roles in determining the overall energy of 
the many-body system in the gas phase. This implies that the energy cost to polarize these 
internal degrees of freedom is small compared to the cohesive energy of the liquid phase. 
Since there is no excitonic correlation in the electron-hole liquid phase, the length scale set by 
the Bohr radius does not have a special meaning to the liquid. As long as the oxide thickness 
is comparable to the inter-electron distance, the liquid phase remains intact except that its 
cohesive and surface energies are reduced. The liquid can be viewed as an entanglement of 
two normal Fermi liquids of electrons and holes |^. Its internal degrees of freedom couple 
strongly to the dynamics of the system. It is a singlet of the spin and valley degrees of 
freedom. In a disordered system, the liquid phase further lowers its energy by occupying the 
valleys of the disorder landscape. It percolates through the sample in a fashion determined 
by the competition among its cohesive energy, its surface energy, and the disorder potential. 
This gives rise to a metal-insulator transition at the percolation threshold in the classical 
limit at zero temperature. Apparently, quantum tunneling between separated liquid regions 
or a finite temperature destroys this transition. This is illustrated in Fig. 0. Along the 
n-axis, i.e. at zero temperature and when the coupling constant determining the phase 
breaking mechanism is infinity, there is a second order phase transition at the percolation 
threshold n = n^. Away from the n-axis, the transition is destroyed by a finite temperature 
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or quantum tunneling. However, in the region of temperature and coupling constant close 
to the critical point, the physics of the system is influenced by the very existence of the 
critical point, similar to the situation in a quantum phase transition [|1^,|11[]. Therefore it is 



the crossover behavior near the critical point that is observed in the experiments. 

We now discuss the related experiments in Si MOSFETs and show how they are consis- 
tent with the consequences of our proposed theoretical framework. 

Temperature dependence of the resistivity — One of the most striking features of the exper- 
iments is that the temperature dependence of the resistivity changes from insulator-like 
to metal-like as the average electron density increases. This is readily understood in our 
theory. At temperatures below the critical temperature and at a low enough average den- 
sity, the system is in the two-phase regime. When the average density is below the threshold 
density of the percolation, the conduction in the system is dominated by the phonon-assisted 
hopping through the localized gas phase. At low enough temperatures, this gives the famous 
Coulomb-gap behavior |]I2| 



P 



;v^, (1) 



where Eq is an energy scale determined by the localization length and the Coulomb inter- 
action. On the other hand, when the average density is above the threshold density, the 
metallic liquid phase percolates through the whole sample. The sample is able to conduct 
electric current. Using a simple two- fluid model for the conduction, we have 

o- = fiCTi + fgag and fi + fg = 1, (2) 

where // and fg are respectively the fractions of the electrons in the liquid and the gas 
phases, a, ai, and ag are respectively the total conductivity, the conductivity of the liquid 
and the conductivity of the localized gas. At temperatures much lower than the cohesive 
energy Ac of the liquid, we estimate fg to be 

fg^Ae-^. (3) 
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Combining the above equations, we obtain the low temperature behavior of the resistivity 
on the metaUic side 

P = Po + Pie"^, (4) 
where po = l/<^« and pi = Apo(l — ag/ai). This temperature dependence of p on the 



metalhc side was first suggested by Pudalov [|T3l based on completely different interpretation 



of the experiments. It was implicit in the experimental data of Kravchenko et al. 0. It 
was explicitly observed by Hanein et al. [H] in recent experiments on GaAs/AlGaAs hole 



samples with a conducting backgate. For T > Ac, the liquid phase ceases to exist. The 
resistivity of the system depends weakly on temperature. 

Scaling and duality — In a realistic sample, the coupling constant of the phase breaking 
mechanism is fixed by the property of the sample. We need to consider what happens in a 
constant-a plane in Fig. |^. There is a diverging length scale ^ and a vanishing energy scale 



e ~ e as the system approaches the critical point Uc- Using the usual arguments [|15 
in the region close to the critical point, we write the finite temperature resistivity of the 
system in the scaling form 

p(T,n)=p(5/T^/'^), (5) 

where 5 = {n — Uc) and u is the correlation length exponent. The value of u predicted by the 
percolation theory is u = 4/3. This is within the range of the measured values of 1.2 ~ 1.6. 

The duality observed in the experiments can be readily understood by the following 
argument. Because there is no critical point on the constant-a plane at a finite temperature, 
the function p{x) is an analytic function of its variable. Taylor expanding it at a; = 0, we 
have 

p(T, n, + S)= p(°) + p(^)^ + lp^'\^r + ■ ■ (6) 
where p^^^ = p(0), p^^^ = ^^^^ \x=o, and p^"^^ = '^^^2^ \x=o- Immediately we have 

p(T, + 6)p{T, n,-6) = {p^'^f + 0(5^). (7) 



Obviously, this operationally defined duality relation is a generic feature of the critical 
region. It should be interesting to test this by looking at the temperature dependence of the 
coefficient of the 5^ term. 

Electric-field dependence of the resistivity — There are two important features associated 
with the electric-field experiments First, the nonlinear I-V occurs at very small current 
(or electric field). Second, the nonlinear resistivity exhibits scaling as a function of electric 
field in both the metallic and insulating sides. There are two known mechanisms through 



which nonlinearity can occur in a homogeneous system [|l5l. The first one requires that 
the energy scale determined by the electric field to be the dominant energy scale in the 
system. In the critical region, this energy should be larger than the temperature. However, 
in a typical nonlinear I-V experiment in Si MOSFETs, the nonlinearity occurs at electric 
fields as low as 0.25 mV/cm 0. Using a phase breaking length of lOOOOA, we estimate the 
nonlinear electric field energy scale to be about 0.3mK, much smaller than the temperature 
220mK. No nonlinearity should result from this mechanism. The second mechanism is due 
to Joule heating. Joule heating raises the effective temperature of the electrons according 
to Tg ~ E'^/^ |TB|. The experiments in Ref. ||^ indicate that Tg is about 2K at an electric 
field of 50mV/cm. Using an electric field of 0.25mV/cm, we estimate Tg to be 150mK. This 
is lower than the lattice temperature 220mK. Therefore Joule heating is not significant at 
such an electric field and no nonlinearity should result. The above arguments leads us to 
believe that the nonlinear I-V is caused by the intrinsic inhomogeneity in the system. 

The scaling behavior appeared at large electric fields can be readily understood following 
the arguments in Ref. [|15|. Close to the critical point, the vanishing energy scale e ~ e^/^ 



is cut off by the electric field through e ~ eE^. Thus, the zero temperature resistivity can 
be written in the scaling form 

p{E,n)=p{S/E'/'^(^^'y), (8) 

where z = 1 because the vanishing energy scale at the critical point is given by e ~ e^/^. 
Using the classic value of = 4/3, we find ^{z + 1) = 8/3 ~ 2.66, which is very close to the 



experimental value of 2.70 §]. 

Effects of an in-plane magnetic field — When the Zeeman energy of the magnetic field 
is larger than the Fermi energy of the gas phase, it fully polarizes the electrons in the gas 
phase. The portion of electrons in the gas phase at low temperatures can be estimated to 
be 

/. ~ e ^ (9) 

when gfisSzH^ < Ac. Using the two-fluid model above, we find the in-plane magnetic field 
dependence of the resistivity on the metallic side at low temperatures to be 

P ~ Po + Pie^^^^^-^ii/^ (10) 

When the Zeeman energy goes beyond A^, the liquid phase ceases to exist because the 
critical temperature is zero. The system is a fully polarized gas and it is not affected 
by further increases of H\\. Fig. |^ shows this generalized situation of Eqn.(|TO|) for a given 
temperature T. The resistivity increases exponentially with H\\ from its zero-field value 
Po(r) to the saturated value Pm{T) at the critical field H\\^c given by qhbSzHw^c = Ac, i.e. 

,^p^T,H,) = I + 4(lnp™(T) -lnpo(T)) if i/,, < H,,^ ^^^^ 

" I In p„(T) if >i7|l,c 

When T <ti guBSzH^ c, piT,H\\) should be exponential in H\\/T for intermediate values of 
H\\. Using this we have 

lnp„(T)~lnp„ + ^^^^|^. (12) 



The preliminary experimental data is consistent with this prediction. It should be 
interesting to test it in more details. 

Next we look at how the low temperature behaviors of the resistivity on the metallic side 
evolve for different in-plane magnetic fields. Using Eqn.(|TT]) and Eqn.(|T^), we have 

lnp(T, i^ii) = (1 - 1^) Inpo(r) + ^^^^ + i Inp^. (13) 



This equation predicts different behaviors of p{T,H\\), shown in the inset of Fig. ^, for 
< //|| and > H^\ Q, where H^\ q = piH^^ ^/iPo + 2pi). Using the experimental data of 
pi/ Pq — 10 and ify ^ — 20 kOe in the experiments of Simonian et al. ||T^, we predict ify^o — 10 
kOe. This is consistent with the data shown in Figure 4 in Ref. ||T^. For H\\ <^ -f^||,0; the 
temperature at which p(T, H\\ ) reaches a minimum is given by 

This is consistent with the prehminary experimental data It should be interesting to 
test this prediction quantitatively. 

Further experiments — (1) One of the most interesting experiments is to probe directly 
the inhomogeneity in the electron density. This can be achieved through scanning techniques 
using a single-electron-transistor as the probe for the electric field generated by the electrons 
in the two-dimensional electron system. According to our theory, we expect the system to 
phase separate into regions of different densities at low enough average electron densities. 
The regions of the sample occupied by the liquid phase should be easily observable in one of 
the implementations of such scanning techniques |]18[. (2) In our theory, the spin and valley 
indices play similar roles in determining the physics of the system. We expect a strain field, 
which splits the valley degeneracy of the two-dimensional electron system in Si MOSFETs, 
to have effects similar to those of an in-plane magnetic field. This should distinguish our 
theory from those in which the spin degrees of freedom play a non-trivial role in producing 
the observed experiments. (3) Light scattering experiments in such systems can also provide 
useful information. Because of the occurrence of macroscopic inhomogeneity in the system 
below the critical temperature Tc, we should see an increase in the wave-vector-breaking 
components of the light scattering as the temperature gets below T^. 

In summary, we have argued that there is a new liquid phase in the two-dimensional 
electron system in Si MOSFETs at low enough average electron densities. The percolation 
transition of this liquid through the disorder landscape in the classical limit in a disordered 
system causes the metal-insulator crossover behaviors observed in recent experiments. We 
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have made quantitative predictions and proposed new experiments to further investigate the 
physical properties of the system. 

Acknowledgments: We thank S. V. Kravchenko, D. Simonian, D. C. Tsui, Senthil Todadri, 
Steve Simmon, and X. R. Wang for stimulating discussion. We also thank the authors of 
Ref. [|r^ for doing the community a big favor by writing such a beautiful article. X. C. Xie 
is supported by the NSF-EPSCOR program. 



9 



REFERENCES 

[1] E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Phys. Rev. 
Lett. 42, 673(1979). 

[2] S. V. Kravchenko, G. V. Kravchenko, J. E. Furneaux, V. M. Pudalov, and M. D'lorio, 
Phys. Rev. 350,8039(1994); S. V. Kravchenko, W. E. Mason, G. E. Bowker, 
J. E. Furneaux, V. M. Pudalov, and M. D'lorio, Phys. Rev. B51, 7038(1995). 

[3] D. Simonian, S. V. Kravchenko, and M. P. Sarachik, Phys. Rev. B55, 13421(1997). 

[4] D. Popovic, A. B. Fowler, and S. Washburn, Phys. Rev. Lett. 79, 1543(1997). 

[5] V. Dobrosavljevic, E. Abrahams, E. Miranda, and S. Ghakravarty, Phys. Rev. Lett. 79, 
455(1997). 

[6] S. V. Kravchenko, D. Simonian, M. P. Sarachik, W. Mason, and J.E. Furneaux, Phys. 
Rev. Lett. 77, 4938(1996). 

[7] D. G. Tsui, private communications. 

[8] T. M. Rice, in Solid State Physics, Vol.32, 1977. edited by H. Ehrenreich, F. Seitz, and 
D. TurnbuU (Academic Press) p.l, and references there in. 

[9] J. G. Hensel, T. G. Phillips, and G. A. Thomas, in Solid State Physics, Vol.32, 1977. 
edited by H. Ehrenreich, F. Seitz, and D. TurnbuU (Academic Press) p. 87, and references 
there in. 

[10] A. J. Millis, Phys. Rev. B48, 7183(1993); S. Sachdev, in proceeding of the 19th lUPAP 
International Conference on Statistical Physics, Xiaman, Ghina, 1995, edited by Hao 
Baihn (World Scientific, Singapore) P.289. 

[11] S. Todadri, Ph. D Thesis, Yale University, 1997. 

[12] A. L. Efros and B. I. Shklovskii, J. Phys. G8, L49(1975). 



10 



[13] V. M. Pudalov, cond-mat preprint #9707076. 

[14] Y. Hanein, U. Meirav, D. Shahar, C. C. Li, D. C. Tsui, and Hadas Shtrikman,cond-mat 
preprint #9709184. 

[15] S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, Rev. Mod. Phys. 69, 315(1997). 

[16] E. Chow, H.P. Wei, S.M. Girvin, and M. Shayegan, Phys. Rev. Lett. 77, 1143(1996). 

[17] D. Simonian, S. V. Kravchenko, M. P. Sarachik, and V. M. Pudalov, Phys. Rev. Lett. 
79 , 2304(1997). 

[18] Amir Yacoby, private communications. 



11 



FIGURES 



FIG. 1. Phase diagram of an electron-hole system. 

FIG. 2. Effects of a finite temperature and quantum tunneling on the percolation 
transition. 

FIG. 3. In-plane magnetic field dependence of the resistivity. The Inset shows the temperature 
dependence of p for various in-plane magnetic fields: (a) ffy <C -ff||,0) (b) < -^[1,0) (c) H\\ > -ff||,o- 
The temperature is measured in unit of g^iBSzHu c- 
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